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We study the establishment of vortex entanglement in remote Bose Einstein condensates. We 
consider a two-mode photonic resource entangled in its orbital angular momentum (OAM) degree 
of freedom and, by exploiting the process of light-to-BEC OAM transfer, demonstrate that such 
entanglement can be efficiently passed to the matter-like systems. Our proposal thus represents a 
building block for novel dissipation-free and long-memory communication channels based on OAM. 
We discuss issues of practical realizability, stressing the feasibility of our scheme and present an 
operative technique for the indirect inference of the set vortex entanglement. 
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I. INTRODUCTION 

The rich variety of coherently exploitable degrees of 
freedom with which a photonic system is endowed has 
been extensively used in recent years in order to demon- 
strate the building blocks of quantum technology pro- 
tocols including quantum cryptography [1], quantum re- 
peaters [5] , teleportation and quantum computing |3] . In 
this context, the exploitation of orbital angular momen- 
tum (OAM) carried by light is settling as a new and ex- 
citing opportunity for coherent manipulation at the clas- 
sical and quantum level [4 . High density data transmis- 
sion [5] , activation of micro-machines and optical tweez- 
ers [6] are among the most prominent applications of op- 
tical OAM so far. 

In addition, the field of quantum information pro- 
cessing has now started exploiting the additional oppor- 
tunities offered by this photonic degree of freedom for 
communication and manipulation purposes. It has been 
shown that it is possible to create OAM-entangled pho- 
tons by means of a routinely used setup such as sponta- 
neous parametric down conversion (SPDC) [7 . This has 
triggered a plethora of studies on how to generate, ma- 
nipulate and detect non-classical states of OAM [5], cul- 
minating in the demonstration of Bell's inequality viola- 
tion by OAM-entangled two-photon states [5], the intro- 
duction of so-called hyper-entangled states [TU], the de- 
sign of quantum cryptographic schemes based on higher- 
dimensional systems [TT] and the transfer of OAM states 
from light to matter- like systems |12l . 

In particular, the latter scenario holds the potential for 
the realization of experimentally feasible long-time quan- 
tum memories embodied by superfluid rotational states 
of Bose-Einstein condensates (BECs) |T^-[Tnj. The spa- 
tial coherence intrinsic in a BEG allows for a superfluid 
vortex state in which the bosons in the condensate have a 
well defined and quantized OAM, which offers a perfect 
match with rotating photon carriers. Along the semi- 
nal lines traced by the experiments in Refs. |12j . a few 
theoretical proposals for the light-to-vortex state trans- 
fer have been presented [T5HTS] . Here we close the circle 



of these proposals and show that it is possible to cre- 
ate entanglement between two spatially separated BECs 
by transferring OAM from entangled photon resources 
to the condensates. We propose a simple and efficient 
scheme to achieve this goal using experimentally achiev- 
able parameters and routinely produced OAM-entangled 
light resources. On a different level, our study proposes 
a scheme that is able to transfer (with in principle 100% 
efficiency) higher-dimensional entanglement between two 
independent system by means of bilocal interactions, thus 
contributing to an area that is witnessing theoretical and 
experimental interest (see Choi et al. in 2J and Ref. 17 ). 

The paper is organized as follows. In Sec. [H] we intro- 
duce the idea behind our proposal and address the Hamil- 
tonian for the OAM entanglement transfer. The key part 
will be the introduction of an adiabatic Hamiltonian in- 
volving Raman processes where photonic OAM quanta 
are used in order to put a single BEC into a vortex state 
with non-zero winding number. Sec. |III| represents the 
core of our proposal by providing a clear physical picture 
of the process discussed above and assessing the quan- 
titative performance of the light-to-BEC entanglement 
transfer. By starting from an OAM-carrying photonic 
state, we show that significant vortex-vortex entangle- 
ment is established by our scheme for experimentally re- 
alistic parameters. Sec. |IV| describes realistic techniques 
for the inference of such quantum correlations based on 
the reverse of the map addressed here. Finally, in Sec. |V] 
we summarize our findings and open up perspectives for 
further development. A few technical steps which com- 
plement the quantitative analysis presented in the body 
of the paper are contained in Appendix A. 

II. THE MODEL 

Let us start by presenting the model used in order to 
describe the light-to-BEC transfer of entanglement. We 
shall see that a key point in this mechanism resides in 
the collective coupling of the atoms belonging to one of 
the BECs to the respective light field. Together with 
the indistinguishability of the resource photons, this per- 
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FIG. 1. (Color online). Sketch of the proposed setup. An 
OAM-entangied two-photon state is produced by spontaneus 
parametric down conversion (SPDC) of a Gaussian pump. 
Each output mode interacts with a respective trapped Bose- 
Einstein condensate (BEC), which is also pumped by an in- 
tense field with no OAM. The local matter-light interaction 
transfers the OAM entanglement from the field modes to the 
condensates rotational degree of freedom. 

mits us to entangle the two BECs. We consider two 
spatially separated and trapped BECs, each with Nq 
(/ = A, B) ^''Rb atoms and let each of them interact 
with one of the field modes of an OAM-entangled two- 
photon state (see the sketch in Fig. [ij. Such a photonic 
resource can be produced, for instance, by type-I para- 
metric down conversion of a Gaussian laser beam which 
is an OAM-preserving process: the sum of the OAM car- 
ried by the entangled signal and idler mode produced 
by a laser-pumped non-linear crystal equals the OAM 
initially carried by the pump |7j. In this paper we shall 
consider a two-photon state produced by SPDC of a laser 
beam with no-OAM [i.e. prepared in a Gaussian spatial 
mode). This implies that the output modes, here labeled 
as a and /?, carry opposite OAM and enter the state 

|$)„^~Co|lo>a|lo>/3+Ci._i|li)„|l_l)^+C_i.i|l-l)a|ll)/3. 

(1) 

Here, \nk)a indicates an n-photon state populating mode 
a and carrying OAM equal to hk. Moreover, the presence 
of |1o)q|1o)^j in \'^)ai} accounts for a residual component 
withnoOAM and |Cop + |Ci -ip |C_i,i p ~ 1 (having 
neglected components with higher OAM, which are very 
weakly populated at moderate pump intensities). The 
main idea behind our proposal is that the arrangement 
of a locally-assisted OAM-transfer from a light mode to 
the respective BEC would also allow for the transfer of 
quantum correlations, therefore constructing an effective 
entangled channel involving remote matter systems. 

The basic building block for the transfer is an off- 
resonant double Raman scattering process. We con- 
sider each individual atom as a six-level system, shown 
in Fig. [2j The energy scheme comprises a ground-state 




FIG. 2. (Color online). Six-level configuration for OAM 
transfer. We show a schematic representation of the relevant 
energy levels of a single *^Rb atom interacting with Laguerre- 
Gauss (LG) driving fields and classical Gaussian pumps with 
frequency LOp. The ground-state triplet comprises states hav- 
ing angular momentum and ±hl. The excited-state triplet 
is adiabatically eliminated from the dynamics by a double off- 
resonant Raman transition with the one-photon detunings A 
and Aq. Two-photon detunings 5±i are also introduced for 
the stabilization of the entanglement-transfer process. The 
component in Eq. [l] carrying zero OAM [being in a Hermite- 
Gauss (HG) spatial mode] drives off-resonantly the |0) o \E) 
transition. 

triplet made out of a non-rotating state |0) and two other 
states, indicated as | ± 1), having angular momentum 
±h. The elements of excited-state triplet |e) and |e ) are 
linked to | ± 1) by two classical pumps, while the field 
modes in \^)ap drive the |0) -H> e, e') transitions. In 
what follows, A and Aq indicate the one-photon Raman 
detunings which are set by appropriate chosing the fre- 
quencies of the driving fields. The classical pumps are 
taken to have a Gaussian spatial profile so that photons 
scattered in the |e) O |1) and |e') -;-> | — 1) transitions 
carry no OAM. Together with the conditions on the OAM 
properties of Eq. ([T]), this ensures that an atom under- 
going the two-photon Raman transition from state |0) to 
I ± 1) (as shown in Fig. [2]) acquires an OAM exactly equal 
to ±h. 

We now introduce the second-quantized matter field 
operators ipij{r) obeying the bosonic commutation rules 
[i^i4r),i;l^{r')] = Si,jS,.jS{r - r'). Here I,J = A,B 

are labels for the BECs while i,j — Q,±l, E,e,e refer 
to the atomic states. As for the photonic part of our 
system, standard commutation relations [c„. k,c^ i.'] — 
^ni,nj^k k' involving the creation (annihilation) operator 
"^1 fe ( ^ni,k) hold. Here ua = a (ns — (3) refers to the 
photonic mode a (/3) that interacts with condensate A 
(B). Finally, k,k' = ±1,0 refers to the OAM degree of 
freedom of the photon. 

Besides the term describing the energy of the free pho- 
tonic fields, the Hamiltonian of the system consists of the 
following four terms 

H = Ha + Haa + Had + Hap, (2) 

which we now address in detail. The first two terms de- 
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scribe the properties of the trapped BECs and are given 
by 



V, 



(3) 



dr V^ljW V-J y (r) V-zj-'W V-ZjCr). 

(4) 

In all the above, V/ is the quantization volume for the 
light-matter interaction involving condensate /. The 
Vij{r)s are the state-dependent atomic trapping poten- 
tials, 77^ ./ = (Anh^ lm)a^. ., accounts for the coUisional 

energy between two atoms (of mass m) in states j and j 
and a\ ^1 is the corresponding s-wave scattering length. 

As it will be clarified later on, the excited triplet \E, e, e'} 
can be adiabatically eliminated from the dynamics of the 
atomic system. Therefore, by assuming no initial popu- 
lation of these states, we can simplify our treatment and 
take a\ TP — a\ ^ = a^. , = for 7 = 0, ±L The third term 
in Eq. ([2]) describes the interaction between BEC / and 
the quantized field mode n/ and can be written as 



'H-ad= V [xi.QCnj.o dr-0| g(r)V'/.o(r)^„,.o(r) 

I=A,B V JVj 



JVi 



+Xi~iCni,^i J rfr ?/'|_g,(r)V'/,o(r)M,-;(r)j + /i.e. 

(5) 

The coefficients xi,k {k = 0, ±/) are the effective 
dipole moments associated with the transitions depicted 
in Fig. [2] The functions £4/.fc(r) describe the spatial 
shape of the states entering |$)q^ and we choose them 
to be 



:(r)=^ 



e , 



(6) 



coupling between the classical pumps and the BECs 



Hap 



I=A.B ^ 



drip] g(r)'07,/(r)e 



-hQje 



—iujit 



Vi 



_,zk/ -r 



(7) 

The coefficients ftj are the Rabi frequencies for the 
matter-pump interactions. Moreover, we initially pre- 
pare each condensate in the atomic state |0) by means of 
optical pumping techniques, for instance, in such a way 
that the excited triplet can be considered as empty. We 
now proceed with the adiabatic elimination of the ex- 
cited states under the assumption of large single-photon 
detunings A and Aq (with respect to the typical coupling 
rates entering Haa, Had and Hap)- We take the time evo- 
lution due to the fields as faster than the center of mass 
motion of the atom when in one of the excited levels, 
so that we can neglect the free atomic Hamiltonian for 
these states. We move to a proper rotating frame where 
we redefine the excited-state (ground-state) matter field 



operators as i/'/.e — ^i,e& , "0/ 



and 



ipi.Ee' 



(^ 



I,k 



-iuJit 



ipi.k) and the photonic 



operators as ci^k[t) = e.'^^'^ci^k (with k — 0, ±/). Fol- 
lowing the works by Marzlin et al. [13] and Kapale and 
Dowling |15j . we explicitly allow for two-photon Raman 
detunings 5±i(t) = ui — ujp{t) — Lj±i (see Fig. [2]), which 
help in the stabilization of the transfer process(see also 
Ref. [H]). Here huj±i are the actual energies of the rotat- 
ing atomic states. In fact, one can intuitively understand 
the necessity for a time-dependent two-photon detuning 
as a result of the adiabatic elimination of the excited 
triplet and the existence of inter-atomic collisions, which 
change the energies of the atom in time. Therefore, in 
order to achieve efficient transfer of 0AM entanglement, 
we need a "chirped" frequency of the pump fields that 
allows to track and compensate the change of the energy 
levels. The technical details behind the adiabatic elimi- 
nation sketched here are presented in Appendix A. 
We thus arrive at the effective Hamiltonian Heff — 

Ha + Haa + Hint for the description of the adiabatic in- 
teraction between light and BECs, where 



1,3 



dr 4,1^ {r)[-^Vi + y,,, (r) + e, {t) ) (r) 



Vi 



2m 



(8) 



here eo is the vacuum permeability. While the first of 
these equations refers to a field mode having a Gaus- 
sian spatial profile, the second describes OAM-carrying 
Laguerre-Gauss beams. We assume that the beam-waist 
W is larger than any linear dimension of the BECs so that 
the Gaussian part of the function reduces to a constant. 
This approximation also ensures the collective nature of 
the interaction between the fields and the atoms belong- 
ing to a given BEC. The last term in H describes the 
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I.k 



Cni.k / dr'^ni.ki'^m fc('^)V'/.o(i") + h.c. 

Vi 



+<,feC„,,fe dr — V),oW'/'/,o(r) 

+ }^Cnj,oCnjfi dr ^ Vi,oW'/'/,o(r)- 



A/,0 



(9) 
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In Eq. ([8]) J = 0, ±1 should be taken, while in Eq. ^ it 
is A: = ±1. Moreover, we have introduced the coupling 
coefficient '^mA'^) = {hn*jXi,k/^i)-(^ni,ki^)e''-^' '' and 
the energies ej(t) such that eo(^) — and e±;(i) = 
h{uj±i — S±i(t) — uj±i). The term describing intcr-particle 
collisions Haa remains identical to Eq. The inter- 

pretation of the form taken by H^ff is straightforward. 
While Eq. ([S]) describes the energy of non-interacting 
matter part of the system, modified by the introduction 
of (5±-related terms, Eq. ^ accounts for the light-matter 
interaction and includes the dynamical a.c. Stark shift ef- 
fect arising from the adiabatic elimination. In particular, 

the first term in Hint describes a three-mode interaction 
where photonic excitations are used in order to perform 
an atomic transition between ground-triplet states. This 
is the key to our analysis on 0AM entanglement transfer 
and the starting point of our quantitative study. 



III. THREE-MODE EXPANSION AND 
LIGHT-INDUCED TRANSFER OF OAM 
ENTANGLEMENT 

A. Bosonic-mode expansion 

In order to fix the ideas and discuss an experimentally 
relevant case, we consider the anisotropic harmonic trap 
potential Vi{r,z) = (m/2)(a;^r^ +uj'^z'^) for each of the 
BECs used in our proposal. Here uj^ (w^) is the frequency 
of the trap along the longitudinal (radial) direction. In 
order to be able to neglect any longitudinal excitations, 
we assume uj^ 3> Wr, so that each BEG is confined in 
a pancake-like structure. In the limit where the inter- 
atomic collisions are very small, the cylindrical symmetry 
of the problem allows us to describe the centre-of-mass of 
one atom in a BEG by means of the set of eigenstates {9 
is the angular coordinate of a cylindrical reference frame) 
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TT 4 ar\/a. 



(10) 



with associated eigenvalues Ei — hwz/'2, + h{\l\ + l)cjr. 
Here the characteristic lengths of the harmonic 

motion along the longitudinal and radial direction. The 
description provided by the eigenfunctions ( 10 1 remains 



valid under the assumption of dilute BEGs, so that their 
ground states result from the simple tensor product of 
the single-particle states (/)/ o(?', 2). This is a good ap- 
proximation as long as N^a{ ^, <^ az [IH]j implying that 
the scattered part of the single-particle wavefunction con- 
tributes with only a small correction to the wave- function 
of the non-interacting case. In order to provide a bet- 
ter picture of the anticipated three-mode interaction de- 
picted in Eq. ([9]), we now define new bosonic opera- 
tors fc/j and bj ■ for the matter-like part of our sys- 



tem as (omitting the tilde from now on for readability) 
?/'7j(r) = (j)ij{r)bij. By using the orthogonality of the 
(j>ij{r) and the commutation relations valid for ip{r)^s, it 
is straightforward to find that [67. i, b\ ^] — dij. Inserting 
these definitions into the effective Hamiltonian we obtain 
a much simplified and self-evident picture of the process 
through 

Ha = ^{Elfi + \k\hUr + €k{t))b\J)I^k, 



I,k 

Hint= ^ {9l.kb\.k^I..0Cni,k + h.c.) 

I,k=l-l 

+ 'b\fibi,o^Pi,kcijA^i,k, 

I,k 



(11) 



where each coefficient can be expressed in terms of the 
non- interacting wave-functions as shown in Table |l] 
The effect of the light-matter coupling is now manifest: 
besides the a.c. Stark shifts proportional to pjj, Hint 
consists of a scattering process at a rate gi^k where the 
annihilation (creation) of a photon of angular momentum 
k is accompanied by the Raman transition |0)/ — \k)j 
{\k)i — ^ |0)/). Such a mechanism, which would deter- 
mine a perfect transfer of OAM from the light resource 
to the BEGs, is disturbed by the inter-atomic collisions in 
Haa and should also take into account the modifications 
induced by the dynamical shifts in Ha,int- The creation 
of inter-BEG OAM entanglement is thus a trade-off be- 
tween these various processes. The task of the next Sub- 
section is precisely the quantitative assessment of such a 
trade-off. It is worth remarking here that, in virtue of the 
definitions of efc(t) the term |fc|/iWr + efe(i) appearing in 
the energy of the rotating states with j = ±1 is explicitly 
dependent on the two-photon Raman detunings 6k (t) and 
takes the form h{ujk — Sk —<^k + l^l^^r)- When the interac- 
tions considered in our scheme are included, the energy 
levels are shifted so that the shift — |A:|aJr 7^ is in 
general non-zero and, possibly, time-dependent. In what 
follows, we shall assume that such shift occurs linearly in 
time. 



TABLE I. Coupling rates in the effective Hamiltonian after 
the introduction of the effective matter-like bosonic opera- 



tors [see Eqs. (Ill] and their expressions in terms of the non- 



interacting atomic wave-function for a pancake-like potential. 



Coefficients 



Corresponding expression 



£x,o /^^ dr0|,o(r)(-^V? + V7,, (r))0,,„(r) 

Pi,o Jv, dr ^/^^ 1 (j>i,o (r) ] 

Pi,i 



J.J 



Jvr '^'^ a" \<t>ifi{r)\ 

V' .L dr|07,,(r)|2|0,y(r)|2 
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B. Entanglement transfer process 

As discussed in Sec. |TT] we assume an initial prepara- 
tion where the atomic excited triplet is empty and all 
the atoms in each BEC populate |0)/. We indicate such 
a collective atomic state as \Nq)j, which condenses in- 
formation on the population of level |0)/ and | ± 1)/. On 
the other hand, the OAM-entangled photonic resource 
is taken as prepared in the experimentally realistic state 

|$z)a/3 = (l/V3)(|lo, lo)c.p + 111, l-l)a/3 + |1-1, ll)a/3). 

(Here we put the label Z which refers to the state 
obtained by Zeilinger and cooworkers) In Refs. [71 [5] 
it was shown that two-photon multidimensional OAM- 
entangled states generated by means of SPDC can be 
effectively distilled into states very close to |$z)a/3, 
thus making the contributions coming from states hav- 
ing higher 0AM negligible [20]. This effectively makes 
the Hilbert space spanned by photonic 0AM states iso- 
morphic to that of a spin— 1 particle, or qutrit, so that 
|$z)q/3 describes a maximally entangled two-qutrit state. 
The initial state is thus taken to be aboB ~ 



\NQ,NQ)AB\'^z)ai3, whose dynamics under Eqs. (Ill is 



now evaluated in a rotating frame defined according to 
Eq. (A-2). 

It is straightforward to verify that, when starting 
from ABafi as given above, the evolved state 

ABafi obtained upon use of the effective Hamilto- 
nian in Eq. ( 11 ) lies entirely in a nine-dimensional sector 



of the Hilbert space spanned by the states 
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-1,1- 


i)aK-1,1i)s|0)„ 





(12) 



The notation used here is such that |A^g 



s,Sk)i indi- 
cates a state where s atoms populate an atomic eigen- 
state of the angular momentum with eigenvalue kh while 
Nq—s atoms populate the state |0)/ having zero 0AM. It 
is worth stressing that the number and structure of the 
states involved in the evolution of a given initial state 
strongly dependent on the total initial angular momen- 
tum carried by the latter. In fact, our effective Hamilto- 
nian preserves the total light-matter 0AM. This property 
implies that the dynamically evolved state of the system 
should be written as 



mt)) 



ABaP 



(13) 
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FIG. 3. (Color online). Evolution of the state probabilities 
|/i(t)|^ against the interaction time t (in lis) for Nq = 10^ 
atoms per condensate. Here f ^ ./ = 5nm, f2/ = 27.5kHz, 
Xi,±i = 1.18kHz, A = 90kHz and the two-photon detunings 
having the functional form S±i{t) = 257/(1 — Qit/2) — ut- 
The trap frequencies are ujt = 70Hz and = 500Hz. The 
photonic resource is tailored at a wavelength of 702nm at an 
angle of 4° off the initial pumping gaussian beam's axis 



with numerical coefficients fi{t) such that |/i(OP = 
1. The analytic solution of such a dynamics is a 
formidable problem and we thus resort to a numerical 
investigation in order to infer the behavior of /i(t)'s. 
To find the coefficients {/i(i)}, we have numerically 
solved the Schrodinger equation using the Hamiltonian in 
Eqs. ( 11 1 [T5^. We have explored a wide range of param- 
eters, including the case where the system is symmetric 
under the exchange of the two BECs, finding qualita- 
tively similar results. In the following, we concentrate on 
the symmetric case and use the parameters listed in the 
caption of Fig. (|3| , which shows that a complete transfer 
of 0AM from the photonic state to the BECs is possi- 
ble, in analogy with the semiclassical case approached 
in Refs. [T5HT5] . The (dashed) green curve shows the 
temporal dynamics of the probabilities |/i,2(i)P whereas 
the (sohd) yellow ones depict |/7,8(i)P for a given set 
of the relevant physical parameters and a specific choice 
for the functional form of the chirped two-photon detun- 
ings. These two sets of probabilities are almost mutu- 
ally mirror symmetric. Damped oscillations are super- 
imposed to a monotonic behavior induced by the com- 
pensation arising from the chirped detunings (in our case 
5±i{t) = 2f]/(l - Vlit/2) - uJt). The low-lying (dotted) 
red curve is for 1/3,4,5,6^, whose corresponding states 
only marginally contribute to the evolution of the sys- 
tem. Finally, the (dot-dashed) horizontal blue line shows 
the probability |/o(OP, which does not change in time as 
l^'o) is an eigenstate of the effective Hamiltonian. As we 
shall see this gives a lower limit for the population of the 
ground state of the reduced density matrix for the two 
BECs. In fact a plot of the populations of the BECs re- 
duced density matrix, pAB{t) = Tra/3(|5'(i))ABa/3(^'(i)|) 
shows that the state \N^) a\N^) b has always a finite 
occupation probability. This is due to the unavoidable 
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presence of photons carrying no 0AM that continuously 
project the two BECs onto their ground states. 

The (solid) yellow curve in Fig. Ill shows the popula- 
tion of the state \B) ^{\N^-1,1i)\n^ -l,l_i) + \N^'^ 
l,l_i)|A^o^-l,li))/V2 (wc have omitted the BEG la- 
bel as no ambiguity exists) which shows 0AM entan- 
glement between the BECs. Clearly, the 0AM entangle- 
ment transfer generates quite a large component of \B) in 
the reduced two-BEC state. This arises from the 0AM- 
carrying components in the photonic resource and the ef- 
ficiency of the population-transfer process. The introduc- 
tion of this state allows us to draw a clear and compact 
picture of the asymptotic form of the map Mt transform- 
ing photonic 0AM entanglement into matter-like one via 
bi-local far off-resonant double Raman coupling. By ne- 
glecting the very small components associated with the 
remaining excited two BEC states ((dotted) red line in 
Fig. |4]) and collecting the remaining terms into a diagonal 
density matrix, this is approximately given by 



lim Mti\N^.Ni^)AB{No^,No^\) 

^i2\B)AB{B\ + \N,^,N,^)AB{N^,N,^\), 



(14) 



Such a formal asymptotic map also explains that we can 
formally infer the properties of the reduced two-vortex 
density matrix by treating it as the state of two (in gen- 
eral entangled) qutrits. 



C. Assessment of entanglement 

We are now in a position to quantitatively estimate the 
amount of vortex entanglement set between the BECs. In 
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FIG. 4. (Color online). Time behavior of the populations 
of the atomic states in pAB{t) = Trafi{\^{t)) ABafi{'^{t)\). 
The same parameters as in Fig. [3] have been used here. The 
(dashed) green curve is for \Nq ,Nq)ab, the (solid) yellow 
one is for state \B). The low-lying (dotted) red curves repre- 
sent the probability that the remaining two-BEC basis states 
are excited. The incomplete population transfer from \0)i to 
I ± l)i (/ = A, B) is due to the zero-OAM terms in |$z)c/3. 



order to tackle this point, our approach will be twofold. 
First, we study the time evolution of the linearized en- 
tropySL{pAB{t)) = (9/8)[l-Tr(pAB(t))] M of the BEC 
density matrix. As Sl is a good measure of the purity 
of a state (it achieves for perfectly pure states and 1 
for statistical mixtures), this will give us an indication of 
the residual entanglement set between the photonic and 
matter-like part of the system: as the dynamics set by 
Hef f is unitary, the entanglement initially present in the 
photonic state has to be conserved when the whole state 
of the matter-light system is considered. Needless to say, 
such entanglement can be transferred from the photonic 
subsystem to the atomic one and/or vice versa and a 
transient can well exist where the two are almost separa- 
ble. Such a situation would be witnessed by a small value 
of Sl and correspond to either large or small values of 
BEC entanglement, with minimal photons-atoms quan- 
tum correlations. We have determined the form of Sl as 
a function of time, which is shown in Fig.[5j As expected, 
at exactly the time when the populations of the 0AM car- 
rying atomic system become non-zero, the linear entropy 
changes its behavior, signaling a maximum of mixedness 
of the light-matter state. This simply implies that for 
t G [400, 500]/j,s the two subsystems are correlated in a 
nonclassical sense. If time increases further, Sl evidently 
decreases, witnessing a reduction in the light-matter en- 
tanglement. Because of the conservation of entanglement 
discussed above, this is the region we are interested in, 
as it could well be the case that in this long-time win- 
dow significant inter-vortex entanglement is set at the 
expenses of the initial all-optical one and the transient 
matter-light correlations highlighted here. 

We quantitatively confirm our expectations by study- 
ing the negativity \12\ [23], an entanglement measure 
based on the violation of the "positivity of partial trans- 
position" (PPT) criterion for separability of a state. Neg- 




800 



t (//s) 



FIG. 5. (Color online). Linear entropy Sl{pab) against time 
t (units of /is). The same parameters as in Fig. [s] have been 
used here. 
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ing N{pab) = 1] is not reached. We stress that this is 
not a Hmitation of our scheme but, on the contrary, an 
effect of the zero-OAM component in the photonic re- 
source Such a detrimental contribution can be 
removed from the BECs reduced state (and its proper- 
ties) by resorting to an "active" approach where, instead 
of discarding the state of Ught after the interaction with 
the condensates, we properly post-select its state. Upon 
inspection of I^'t.s) in Eqs. (12 1, it is straightforward to 



see that state \B)ab is associated with modes 1 and 2 in 
the vacuum state. On the other hand, the entanglement- 



400 

t {lis) 



spoiling component 



)ab would bring about pho- 



FIG. 6. (Color online). Negativity N{pab) against the in- 
teraction time t (units of ^s). The same parameters as in 
Fig. [3] have been used here. We obtain a maximum of entan- 
glement in correspondence of the range of decreasing trend of 
<S'i(pAs(i)), which witnesses a larger purity of the two- vortex 
state and smaller quantum correlations between light and 
matter. Inset: Negativity against t for Nq = lOiVo' = lO" 
and SIb = 10 Q,a = 27.5 kHz. Other parameters as in Fig. [3] 



"'^Hons in both the modes. It is therefore sufficient to use a 
standard Geiger-like avalanche photo-detector per mode, 
which discriminate the vacuum from the presence of any 
non-zero number of photons in a field, in order to operate 
the optimal post-selection of the BECs state: by register- 
ing no click at both the photo-detectors, we project the 
state of A and B onto the maximally entangled state \B). 
It is in fact worth stressing that, by effectively excluding 
the possibility that the atoms occupy state |0)/, the post- 
selection procedure further reduces the dimension of the 
relevant Hilbert space spanned by each vortex state to 
a bidimensional one, thus leaving us with two effective 
qubits. 



ativity is defined as 



N{pAB{t)) = -2Y,K 



(15) 



where AjT are the negative eigenvalues of the partially 
transposed density matrix with respect to one of the BEC 
systems. The results are shown in Fig. [6] As expected, 
the region of large inter- vortex entanglement corresponds 
to the range of interaction times where the linearized en- 
tropy decreases towards a steady-state value. We point 
out that the wavy behaviour of the curve in Fig. (|6]), 
as well as in the other plots of the paper, is due to the 
inter-atomic scattering. This is confirmed by the the plot 
shown in the inset of the figure, where an asymmetric 
case has been studied. The comparison between sym- 
metric and asymmetric case shows that a mismatched 
number of atoms in the two BECs results in a change 
in the oscillatory behavior of the curve describing the 
time evolution of the transferred entanglement. On the 
other hand, mismatched Rabi frequencies only determine 
a change in the temporal scale of the entanglement dy- 
namics. The degree of transferred entanglement is only 
mildly affected, which demonstrates the robustness of our 
protocol to such effects. The efhciency of our protocol 
does not depend on the assumption of symmetry under 
the exchange of the condensates and is retained in a wide 
range of the relevant parameters. Therefore, for the sake 
of convenience and without affecting the generality of our 
results, in what follows we restrict our attention to the 
symmetric case. 

As discussed above, quite a large value of entangle- 
ment is set between the vortex states of the two con- 
densates, although a maximally entangled state [achiev- 



IV. DETECTION OF VORTEX 
ENTANGLEMENT 

In this Section we describe a method for the detection 
of the vortex entanglement created by the process above. 
Given the low-excitation level of our protocol, the usual 
matter-wave interference is not helpful and we instead 
propose an approach based on the inversion of the pro- 
cess addressed here for light-to-BEC entanglement trans- 
fer. After generation of a two-vortex entangled state (as 
described in Sec. |IIip , the OAM-transferring interaction 
should be stopped. We thus assume that the pump fields 
have been turned off (or set far off-resonant with respect 
to the frequency of the transitions they guide) so as not 
to perturb the entangled states of the vortices. The time- 
reversal nature of our protocol makes it intuitive to un- 
derstand that, if we now reinstate such pumps, photons 
will be scattered into two Laguerre-Gauss modes at the 
frequency of the |0)/ O |e, e') transition, thus writing 
back the two-vortex state onto light fields. One can then 
apply state-property reconstruction techniques, includ- 
ing testing Bell's inequality violation for bipartite states 
of effective three-level particles ^25^. However, the success 
of such tests is usually very sensitive to the form of the 
state under scrutiny and the level of non-ideality affect- 
ing it. We thus resort to a specific and quite promising 
way to infer the properties of the state we have gener- 
ated based on Wigner function reconstruction, which is 
possible by using computer generated holograms [7] and 
homodyne-like measurements |26) . 

Theoretically, the Wigner function for the 0AM state 
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of a photon has been defined in the discrete cyhnder 
Z X Ci (Ci is the unitary circle) representing the phase 
space for the 0AM operator and its canonically conju- 
gate operator 9. In Ref. [21] it has been shown that the 
study of the Wigner function for an 0AM state gives in- 
formation both on the various 0AM eigenstates involved 
in the description of the state and on their relative phase. 
We define the two-mode Wigner function as 

IbiOb) — Tr[w(Z^, (g) w(/_B, (16) 

where we have introduced the kernel w{lj,9j) (/ — 
A, B) mapping quantum states in phase space [21] in 
a way completely analogous to the continuous position- 
momentum phase space. In Fig. [t] a) we compare 
the Wigner function VF(0, 6, 0, 0) associated with pAsit) 
when the entanglement transfer has not occurred [wide 
blue circle] to what is achieved at long-enough t, where 
the map Ait has been implemented [inner (red) butterfly- 
like structure]. The difference due to the coherence es- 
tablished between two-mode orthogonal eigenstates of 
the 0AM operator in the entangled state pab is strik- 
ing. Moreover, from Fig. [t] b) and c), we find that 
W{0, 9a, 0, 9b) exhibits oscillations whose amplitude de- 
pends on the relative phase between the states \B) and 
\D). The curves shown in panels b) and c) are associ- 
ated with the same interaction time as in the circular and 
burtterfly-like structures shown in panel a) , respectively, 
which demonstrates that the oscillation amplitudes de- 
pend on the populations of these states. Remarkably, 
the analysis of W{Ia, 6 a, Ib, ^b) for incoherent superpo- 
sitions of 0AM eigenstates results in a flat distribution. 
Therefore, by experimentally reconstructing the Wigner 




FIG. 7. (Color online). Analysis of the Wigner function for 
pAnit). (a): The blue circle shows a projection the Wigner 
function on the unit circle for t = 0, while the red but- 
terfly structure is associated with t ~ SOO^s, where OAM 
entanglement has been transferred, (b) and (c): We plot 
W{Q, Oa,Q, 9b) at the two instants of time considered for panel 
(a). Notice the dilTerent vertical-axis scales in the two plots. 
The visibility of the fringes of interference in the Wigner func- 
tion is an indication of quantum correlations. 



function, one can determine the inference of entangle- 
ment in the associated OAM state. Any deviation from 
flat distributions typical of incoherent superpositions im- 
plies coherence in the bipartite state, although not en- 
tanglement. Methods based on the inverse Radon trans- 
form [27] could then be used in order to achieve full in- 
formation on the state and, eventually, the entanglement 
set by the transfer mechanism. 



V. CONCLUSIONS AND OUTLOOK 

We have shown that vortex states of spatially re- 
mote, non-isotropically trapped BECs can be entangled 
by means of bilocal OAM-transfer processes and quan- 
tum correlated photonic resources. The amount of vortex 
entanglement set by our scheme can be quite considerable 
and appears to be limited only by the zero-OAM carrying 
component in the photonic resource. While such a bottle- 
neck can be actively bypassed by means of post-selection, 
as described in Sec. IIII] we are currently working on a 
modification of our protocol based on the use of other 
forms of the entangled photon-pair [28j. The difference 
between the two-vortex state achieved by our scheme 
and a classical admixture of OAM eigenstates (without 
coherence and, thus, entanglement) can be revealed by 
a straigthforward state-retrieval process and the recon- 
struction of the OAM-state Wigner function. We believe 
that the superfluid phase of a BEG, together with virtu- 
ally frictionless rotational states of light-induced vortices, 
can be reliably exploited in order to set a promising sce- 
nario for the storage of quantum information and the 
distribution of quantum correlated channels for commu- 
nication. 
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APPENDIX A 

In this Appendix we present the calculations which 
lead to the expression of the effective Hamiltonian in 
Sec. [nj The first step is the adiabatic elimination of the 
excite triplet {|i?),|e),|e )}, which wc perform by defin- 
ing the rotating-picture new matter field operators 4'i.e = 

?A7,ee-'"'*, ipjy = ^/,e'e"*"-'* e V-z^b = ^/.se-'""* with 

V'/,ei V'/ e' } taken as time-independent. Such op- 
erators are then used in the Heisenberg equations for the 



dynamics of the excited-triplet states. One finds 



L „iki-r -i{uji-uj,)t , 



A/ 

hAj 



Cni,i V'/,o(r), 



(A-1) 



where we have defined the single photon detunings 
A/,0 = Wo - w/^B, A/ = - w/^e. In order 
to explicitly include the chirped two-photon Raman 
detunings, which are crucial in the stabilization of 
the transfer process, we define new field operators in 
a rotating frame defined by the Hermitian operator 

O^J2i.k^i (■^U(i",^)'0/,fc(r,O - cl^j^Cn^^kY Explicitly 

^/,fc(r,i) = e^O*V>/,fc(r,f)e-'0* = e-^^^^V'/.felr, t), while 
the photonic operators become Cnj^k = e*'^*c„j-_fce"*'^* = 
Defining the free Hamiltonians Hq''^ — 



-h^W]/{2m) + yf,fc(r) and with the help of Eqs. (|A 
[T]) we finally get 



d = 



+ 



hAi 

Ai 



(A-2) 



(A-3) 



A/ 



Vi 



dr — ■(/'|,o(r)V'/,o(r)Jc„,,o 



hAin 



(A-4) 



rfr V'|,oWV'/,o(r)]cn^,! 



(A-5) 

From these expressions it is straightforward to define the 
effective interaction Hamiltonian H^ff used throughout 
the paper. 
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